Quantum and Classical Orientational Ordering in Solid Hydrogen 



OS 

o> 



(N 

> 

oo 
O 
in 
o 

\o 

c3 



o 
o 



13 



I. I. Mazin, Russell J. Hcmlcy, A. F. Goncharov, Michael Hanfland, and Ho-kwang Mao 
Geophysical Laboratory and Center for High-Pressure Research 
Carnegie Institution of Washington, Washington, D.C. 20005-1305 

We present a unified view of orientational ordering in phases I, II, and III of solid hydrogen. Phases 
II and III are orientationally ordered, but the ordering objects in phase II are angular momenta 
of rotating molecules, whereas in phase III the molecules order themselves. This concept provides 
quantitative explanation for the vibron softening, libron and roton spectra, and the increase of the 
vibron effective charge in phase III, as well as a framework for understanding the topology of the 
phase diagram and ortho-para state at high pressure. The effective charge and the IR and Raman 
vibron frequency shifts are all linear in the order parameter in phase III. 



In the last decade, a wealth of information about the 
phase diagram of hydrogen at high pressures has been 
collected. However, even a qualitative understanding of 
the phase transformations in the solid is not yet in hand. 
Three phases arc known in the experimentally accessible 
range of pressures, which extend to 300 GPa and span 
over an order of magnitude in compression (Fig. |l|). The 
high-temperature phase consists of a closed packed lattice 
of freely rotating and, on average, spherically symmet- 
ric molecules jD. The two low-temperature phases have 
lower symmetry, which implies at least partial ordering 
of the molecules |l]-[|] ■ The properties of the two phases 
are so distinct, however, that more general concepts than 
mere crystallographic dissimilarities must be invoked to 
understand the origin of various phenomena they exhibit. 
The most unusual features include a dramatic increase 
in the infrared activity of the main vibron, a strong vi- 
bron softening, an unconventional geometry of the phase 
boundaries, and drastic changes in rotational excitations 
0-0]. Here we show that these disparate observations 
can be understood in terms of the concept of quantum 
versus classical orientational ordering in the dense solid. 

Hydrogen at low temperature and pressure forms the 
only molecular quantum solid. As such, the excitations of 
the freely rotating molecules in phase I can be described 
by the rotational quantum number J. This contrasts 
with heavier molecular crystals, where molecular rota- 
tion is substantially hindered even at low pressures ||. 
The properties of H2 and D2 crystals greatly depend on 
the ground state of the constituent molecules (J=0 or 1), 
which in turn is determined by their total nuclear spin. 
Here we consider hydrogen solids having a substantial 
fraction of spherically symmetric (J=0) molecules (para- 
H2 or ortho-D2). Solids containing a sufficiently high 
concentration of J=l molecules (e.g., 0-H2 / P-D2) trans- 
form to a cubic phase (space group Pa3>) even at ambient 
pressures (e.g., 3.8 K) M. In this transition, the direction 
of angular momenta of individual molecules orientation- 
ally order, a transformation driven by quadrupolar in- 
teractions, where the ordered phase is the lowest-energy 
configuration for classical quadrupoles || . 



A different kind of orientational ordering occurs in 
phase II. In this case, solids consisting of spherically sym- 
metric molecules (p-H2 or 0-D2) transform to ordered 
(broken symmetry) phases at ~110 GPa in P-H2 j| or at 
~28 GPa in o-D 2 @ at T -> K (Fig.0). This transition 
was interpreted as arising from increasing intermolecu- 
lar interaction that results in mixing of higher angular 
momenta (e.g., J = 2, 4 for p-H.2) into 
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FIG. 1. Phase diagram of hydrogen at megabar pressures 
obtained from Raman and infrared measurements p)fl[12|. 
The open triangle and diamonds are from Ref. Q 

the ground state molecular wave function and impart- 
ing a finite angular moment for the molecules The 
non-zero angular momenta can then order, with molecu- 
lar centers forming a lattice derived from hep Q]. Even 
more exotic ordering schemes may exist for normal H2 or 
D2 0] ■ We will refer to these low-pressure ordered phases 
collectively as phase II. Here, ortho-para distinctions are 
valid in the sense that the wave function can be defined 
for individual molecules and is either even or odd. 

One expects that at still higher pressure intermolec- 
ular interactions become so strong that all higher J 
states have the same weight; that is, the molecules be- 
have classically and can orientational order as classi- 



1 



cal rotors. Such a pressured-driven transformation be- 
tween two kinds of orientationally ordered states would 
be a transition between a quantum and classical crys- 
tal. This would be unique in relation to known quantum 
crystalline transitions (e.g., solid He) in the sense the 
quantum/classical distinction considered here is associ- 
ated with rotational degrees of freedom. As discussed 
above, spectroscopic changes at the II-III transition near 
150 GPa are significantly larger than those due to ori- 
entational ordering at the I-II boundary. Not surpris- 
ingly, the higher pressure transition has been the focus 
of considerable theoretical study recently (e.g., [p|-^0t). 
We show that new and previously reported data, includ- 
ing the orientational order parameter, vibron effective 
charge, changes in roton and libron spectra, geometry of 
the phase diagram, and evolution of the ortho-para state, 
point to the transition to phase III as being such a trans- 
formation. This provides a simple and transparent model 
for orientational ordering in phases I, II, and III. 

1. Order Parameters and Infrared Intensity. We first 
consider the increase in IR vibron absorption associated 
with passage into phases II and III (Fig. ||). At 167 
GPa and 85 K, the effective charge of the vibron is q* ~ 
0.032e Q; it increases with pressure and reaches 0.037e 
at 230 GPa. In phase II, in contrast, q* is < 0.004e at 
140 GPa and the same temperature, whereas phase I is 
characterized 
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FIG. 2. Temperature dependence of IR vibron spectra, (a) 
I-II transition, (b) I-III transition. 



by a broad, weak disorder-induced band jl 12 F|. The 
origin of these differences can be understood from analy- 
sis of the temperature dependence of q* . We find a strik- 
ing correlation between the temperature dependence of 
(q*) 2 (integrated intensity) and Av, the frequency shift 
with respect to its value at the transition point (P C ,T C ) 
for phase III; i.e., (q*) 2 oc Av 2 over the entire P-T range 
investigated (Fig. Bp,). The observation of q* oc Av can 



be interpreted to mean that both quantities are linear 
functions of a scalar order parameter, Av oc r), and q* oc r\ 
pdj . Recent measurements of the Raman vibron in phase 
III (for D 2 ) |H| showed that the temperature dependence 
of the frequency shift (proportional to 77) can be described 
by a Maier-Saupe model, which characterizes the orienta- 
tional ordering of classical rotors and initially was derived 
for liquid crystals. We show here that the same is true 
for the IR intensity and both IR and Raman shifts [jl]j|] 
(Fig.^b). In contrast, the normalized order parameter in 
phase II is qualitatively different, being much steeper as 



a function of T/T c (cf., Ref. JT3J). 

2. Magnitude of the Effective Charge. We now show 
that not only the temperature dependence of n, but also 
the magnitude of q* can be understood within the pro- 
posed framework. Previously, the origin of the vibron 
intensity in phase III was examined from various elec- 
tronic standpoints ■ We note here that a condition 
for vibron IR activity is that the two atoms in a molecule 
are crystallographically 
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FIG. 3. (a) Experimental integrated IR intensity versus fre- 
quency shift Av for phase III. (b) Order parameter from IR 
intensity, and IR and Raman frequency shifts at 176 GPa. 

inequivalent |16| . If the molecules have non-zero 
quadrupole moments, a non-zero electric field is then in- 
duced at the lattice sites, which in turn polarizes the 
molecules and creates an effective dipole moment [ jl7| . 
We have estimated the magnitude of the effect for H2 
by calculating the induced field in several well-studied 
quadrupolar lattices. The total field is proportional to 
the quadrupole moment, E — f3Q/a 4 , where a is the lat- 
tice parameter. The component of the field along the H- 
H bond, Eu, is f3\\Q/a 4 . The values for (5 depend on the 
structure and orientation, but for a large class of struc- 
tures they are of order unity [Q . For a static molecule in 
the zero-pressure solid, Q ~ 0.5 a.u. || and near the II-III 
transition a is ~3.3 bohr. At the equilibrium bond length 
(d = 1.5 bohr) the polarizability a 11 of the H 2 molecule 
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is 6.72 bohr~ 3 (R|, which gives an effective static charge 
q = aE/d w 0.025e. To find the dynamic charge, we sub- 
stitute a by d(da/dd). From Ref. [ p^| , we have qd ~ 1.5g, 
which gives q* » 0.06e for the total effective charge, close 
to the experimental value of 0.03-0.04e 20 

This result should be contrasted with the behavior of 
the vibron in "quantum" phase II. In this phase, the 
narrow absorption peak appears on the shoulder of the 
broad disorder-induced band (Fig. ||) ^M. This sharp 
peak is a symmetry-allowed IR-active vibron and signals 
the onset of orientational ordering of angular momenta. 
The maximum intensity of this vibron (q* » 0.004e) is 
an order of magnitude smaller than that for phase III 
(e.g., at 165 GPa) at the same temperature fllTf. Sev- 
eral factors reduce the effective charge in phase II com- 
pared to that in III: (a) In a pure quantum state in which 
{J,Mj} — {1,1}, the expectation value of Q is reduced 
by a factor of (21 + |^l| with respect to that of a 
static molecule, and (b) the polarizability perpendicu- 
lar to the molecular bonds is smaller than that parallel. 
The smaller polarizability not only diminishes the static 
charge by 40%, but also reduces the dynamic charge by 
a factor of three. As a result, the effective charge is re- 
duced by 2-3.5 from the change in (Q) and 2-2.5 from 
the polarizability change, giving a total reduction by a 
factor of 5-10, in agreement with experiment. 

3. Rotons and Librons. The proposed picture fur- 
ther implies that elementary excitations corresponding 
to the angular degrees of freedom in the system (rotons 
and librons) must be very different in phases II and III. 
Indeed, IR and Raman data reveal a striking change in 
these low-frequency excitations upon crossing the II-III 
phase boundary. Fig.^ shows rotational and librational 
mode frequencies from new and earlier spectroscopic data 
up to 230 GPa. No discontinuity is observed in the broad 
roton bands at the I-II transition for normal H2 (110 GPa 
and 85 K). In contrast, the rotons disappear at the II- 
III transition, and are replaced by new excitations in the 
same energy range. This change can be understood by 
considering the difference in libron spectra for quantum 
and classically ordered phases. For the former, J remains 
a good quantum number and the roton excitations char- 
acteristic of the disordered phase are supplemented by 
Mj — > Mji "libron" excitations of the ordered phase, 
as in the Pa3 structure of ortho-rich H2 [^). Likewise, 
phase II consists of rotating molecules with a similar set 
of excitations. In contrast, the molecules in phase III 
behave as classical objects (as in solid N2), and J is no 
longer a good quantum number. The elementary exci- 
tations associated with rotational degrees of freedom are 
neither J — » J' rotons nor Mj — > Mji transitions but 
are anharmonic classical librons, quantized angular os- 
cillations about the equilibrium orientation in which full 
rotation is strongly hindered. Moreover, the phase III 
librons harden dramatically with pressure and extrapo- 
late to zero about 70 GPa below transition (Fig.||). This 



strong pressure dependence is expected for librons and 
can be contrasted with the behavior 
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FIG. 4. Rotational and librational frequencies for phases 
I, II and III of hydrogen (85 K) [§,§Q. The IR data were 
obtained by subtracting the vibron frequency from IR com- 
bination bands measured. A new Raman band in phase III 
having a strong pressure dependence identical to one of the IR 
modes is shown. The dashed line suggests "latent soft mode" 
behavior of the librons. 

of the higher frequency lattice mode, which shows at 
most a small discontinuity across the II-III transition [0 . 
These observations, together with the continuous and dis- 
continuous changes in the vibrons fl5f| , are indicative of 
a first-order, but non-reconstructive, mechanism for the 
transition. 

4- Geometry of the Phase Diagram. The geometry 
of the phase boundaries near the I-II-III triple point 
(Fig.|l|) also fits the picture of a "quantum-classical" ori- 
entational transition. Near the triple point, now estab- 
lished for H 2 and D 2 (l|,|l2|,^3|Jl|] , the I-II boundary is 
nearly horizontal, dT c /dP c < 0.4 K/GPa, while the II-III 
boundary is nearly vertical up to ~80 K, dT c /dP c > 30 
K/GPa in H2; this effect is even more pronounced in 
D 2 p5|. Such a configuration is highly unusual for a 
crystallographic transition. One can use the Clapeyron 
equation, AV/AS = dT c /dP c , to estimate the volume 
change across each phase transition. Assuming an en- 
tropy change DeltaS « 0.2i?/mole (Ref. g, p. 238), 
we estimate DeltaV at the I-II transition to be 0.0007 
cm 3 /mol. An analogous calculation for the I, II-III tran- 
sition gives about 0.05 cm 3 /mol. Indeed, one should ex- 
pect large changes in intermolecular forces for the lat- 
ter because of modifications in quadrupole-quadrupole 
interactions 0. This effect can be estimated for static 
molecules from the electrostatic energy of quadrupolar 
closed-packed lattices, U — 1.5kQ 2 /a 5 , k « 7 |. Dif- 
ferentiating this with respect to volume, we calculate an 
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"electrostatic" pressure of about 8 GPa, which gives 0.05 
cm 3 /mol at 150 GPa p^j . Note that the quadrupolar 
energy is reduced by a factor of 4/25 for rotating 
J = 1 molecules, which is in qualitative agreement with 
the much smaller volume change at the I-II transition. 

5. Ortho-Para State. The "quantum-classical" transi- 
tion associated with phase III has important implications 
for the evolution of ortho-para states as a function of 
pressure. There is evidence for ortho-para distinguisha- 
bility under moderate pressures, but its persistence to 
the megabar range has been a subject of debate 0. The 
above analysis naturally leads to the conclusion that the 
wave function of the crystal in phase III can no longer 
be factored into individual molecular wave functions, so 
single- molecule ortho-para distinguishability is lost. In 
contrast, despite the mixing of J states in phase II, the 
notion of parity of the wave function of nuclei in an in- 
dividual molecule is valid. This is supported by recent 
measurements for D2: in phases I and II significant dif- 
ferences in IR and Raman bands are observed for samples 
starting out as pure ortho versus normal @; in phase III, 
however, the number and frequencies are identical and 
independent of sample history |l^]L3[ ]. 

In conclusion, we suggest that orientational ordering in 
the higher pressure molecular phase of hydrogen (phase 
III) is qualitatively different from that observed at lower 
pressures (phases I and II) . The molecules in phase III are 
orientationally ordered in the sense that the time aver- 
age of the molecular bond direction is non-zero, whereas 
in phase II the ordering objects are axes of quantization 
of angular momentum. This concept provides a natural 
explanation of the striking IR vibron activity and other 
properties of phase III. Further quantitative treatment 
should take into account detailed structural changes oc- 
curring at these transitions as well as the effects of zero 
point motion (e.g., Ref. ||]). 

We thank M. Li and J. F. Shu for assistance, and N. 
W. Ashcroft and J. H. Eggert for comments. This work 
was supported by the NSF and NASA. 



[1] H. K. Mao and R. J. Hemley, Rev. Mod. Phys. 66, 671 
(1994). 

[2] R. J. Hemley and H. K. Mao, Phys. Rev. Lett. 61, 857 

(1988); ibid, 63, 1393 (1989). 
[3] H. E. Lorenzana, I. F. Silvera, and K. A. Goettel, Phys. 

Rev. Lett. 63, 2080 (1989); ibid, 64, 1939 (1990). 
[4] M. Hanfland, R. J. Hemley and H. K. Mao, Phys. Rev. 

Lett. 70, 3760 (1993). 
[5] J. van Kranendonk, Solid Hydrogen, Plenum, N.Y, 1983. 
[6] I. F. Silvera and R. J. Wijngaarden, Phys. Rev. Lett. 47, 

39 (1981). 

[7] A. F. Goncharov et al., Phys. Rev. B 54, R15590 (1996). 



[8] R. J. Hemley, Z. G. Soos, M. Hanfland and H. K. Mao, 

Nature 369, 384 (1994). 
[9] B. Edwards, N. W. Ashcroft, and T. Lenosky, Europhys. 

Lett. 34, 519 (1996); B. Edwards and N. W. Ashcroft, to 

be published. 

[10] J. S. Tse and D. D. Klug, Nature, 378, 595 (1995). 
[11] R. J. Hemley, I. I. Mazin, A. F. Goncharov, H. K. Mao, 

Europhys. Lett., to be published. 
[12] R. J. Hemley and H. K. Mao, in Elementary Processes 

in Dense Plasmas, edited by S. Ichimaru and S. Ogata 

(Addison-Welsey, Reading, Mass., 1995) p. 269. 
[13] L. Cm, N. H. Chen, and I. F. Silvera, Phys. Rev. B 51, 

14987 (1995). 

[14] We note that the orientational order parameter is not a 
scalar, but a tensor, n a p — (3r a r,3 — <5 a /3), where r is 
the direction of the molecular bond, rj^p has 5 indepen- 
dent components, given by three Euler angles, its magni- 
tude, and its asymmetry. Only if the average orientation 
is temperature independent (or weakly dependent) and 
the asymmetry of the order parameter is negligible, can 
it be reduced to a one-component order parameter, as ob- 
served experimentally. The connection between the scalar 
order parameter, vibron softening, and IR activity is also 
examined in Refs. [B,|ljl. 

[15] A. F. Goncharov et al, Phys. Rev. Lett. 75, 2514 (1995). 

[16] R. Zallen, R. M. Martin and V. Natoli, Phys. Rev. B 49, 
7032 (1994). 

[17] The quadrupole-quadrupole interaction is not the main 
factor driving orientational ordering in phase III (e.g., 
Ref. ||^2|). Nevertheless, as long as molecules are or- 
dered, their quadrupole moments can play an essential 
role in vibron IR activity. 

[18] For example, in the "herringbone" structure proposed by 
Kaxiras et al. [Phys. Rev. Lett. 76, 1138 (1991)], /3's de- 
pend on the tilting angle 9. At 8 — ((3,/3u) = 0, but 
at 9 — 90° (/3,/3y) = (1.3,0.66); in the equilibrium struc- 
ture 9 = 54° and Q3,/3||) = (0.85,0.34). In the Pca2i 
structure p2| , (/3, /3y) depend on the orientational an- 
gles (9,(f>) and vary from zero at 9 = to (1.35,1.35) at 
{9,4>) = (90°, 90°). At the equilibrium orientation §2), 
(6,<j>) « (62°, 50°) (/3,/3||) = (0.3,0); even at T = OK, 
however, zero-point motion forces the average orientation 
to be (6,<j)) w (62°, 90°), where {fi,P\\) = (1.05,0.95). 

[19] W. Kolos and L. Wolniewicz, J. Chem. Phys. 46, 1426, 
(1967). 

[20] The calculated value is expected to be larger because the 
effects of finite temperature and zero-point motion are 
neglected in the model. 

[21] If J > 1, the state with Mj — has the largest absolute 
value of the quadrupolar moment (which is negative) ; for 
J — 1, however, the Mj = state has twice the expecta- 
tion value of Q relative to Mj = ±1 (and is positive). 

[22] I. I. Mazin and R. E. Cohen, Phys. Rev. B 52, R8597 
(1995). 

[23] This is calculated from the slope of the equation of state 
at 150 GPa, which is approximately 160 GPa-cm _3 -mol 
[P. Loubeyre et al., Nature 383, 702, 1996)]. 



4 



